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We calculate the persistent current of ID rings of spinless fermions with short- 



range interactions on a lattice with up to 20 sites, and in the presence of disorder, for 
various band fillings. We find that both disorder and interactions always decrease the 
persistent current by localizing the electrons. Away from half-filling, the interaction 
has a much stronger influence in the presence of disorder than in the pure case. 
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The recent discovery of persistent currents in mesoscopic rings has addressed new in- 
teresting questions on the thermodynamics of these systems. Although such an effect was 
predicted for a long time, the unexpectedly large amplitude of the measured currents lead 
to important interrogations. Among them, the role of e-e interactions is still unclear. It has 
been proposed [I| that the interactions contribute to the average current which is measured 
in a many rings experiments f2|. On the other hand, the importance of the choice of the 
statistical ensemble to calculate average values has also been stressed ||. Although the first 
explanation, based on a perturbative calculation both in interaction and disorder, seems to 
give a quantitative estimate closer to the experiment, it is still too small by one order of 
magnitude and the interaction parameter used in the theory is not well known. 

In addition, for a single ring experiment [|J , the magnitude of the measured current is also 
not understood and, up to now, perturbation theory has failed to explain an enhancement 
of the current ||. It is only when disorder is weaker that experiment and theory seem to 
agree ||, even for non interacting electrons. 

At the moment, the role of the interactions in disordered systems is still unclear and the 
subject is vastly open. It has been recently proposed that, in the presence of interactions, 
the current should be larger than for free electrons, the effect of the interactions being to 
counteract the disorder effect 0. 

The aim of this paper is to describe the interplay between the interactions and the 
disorder on the persistent currents in ID rings. We choose a model of spinless fermions with 
short range interactions on a lattice. Our main result is that a repulsive interaction always 
decreases the amplitude of the current. It is well known that the ID description is certainly 
not the most appropriate one to describe quantitatively experiments which are performed 
in rings with finite width, in the diffusive regime. But our hope is to find numerical results 
which may give indications for a more real situation. 

We describe a chain of ID spinless fermions in the presence of disorder with the following 
Hamiltonian : 

H = -t/2 ^](exp (2m$/N) c\ c i+1 + h.c.) + V ^ ni n i+1 + J2 w i n i C - 1 ) 

i i i 

where W{ are on-site energies and are chosen randomly between -W/2 and W/2 and V is 



the nearest neighbour Coulomb repulsion. In the following we will take t=l, $ is the total 
magnetic flux through the ring (measured in units of flux quantum <3> = h/e) and N is the 
number of sites. (We use the conventional notation for the amplitude of disorder W and for 
the nearest neighbors interaction V. These notations are opposite to those used in ref. 0). 

Let us first recall some physical properties of this hamiltonian without disorder, i.e. 
W=0. In one dimension, for repulsive interaction, a metal-insulator transistion occurs at 
half-filling due to the existence of umklapp processes. However away from half-filling, the 
umklapp processes become irrelevant and the system is expected to be metallic §|. For 
Hamiltonian ( p.l|) , the metal-insulator (Mott) transition occurs at V=l. The system is 
insulator for V > 1 and metallic for < V < 1 . 

This transition can also be described in the spin picture ||. In the case W = 0, the 
model is integrable and is formally equivalent to an anisotropic spin model, as obtained by 
a standard Wigner- Jordan transformation. By this way the new Hamiltonian reads: 

Hxxz = -t/2 £(exp(2m$/iV) Sr +1 + h.c) + V[N/A + £ Sf 5f +1 ] (0.2) 

i i 

For the XXZ model, V=l corresponds to a transition from an XY-model (V < 1) to an 
Ising model (V > 1). V = 1 corresponds to an isotropic Heisenberg system. Note that 
a spin gap opens up for V > 1 . It corresponds to the gap in the charge excitations for 
Hamiltonian (p.l|) characteristic of the insulator. With disorder, W ^ 0, we must add in 
( p.2| ), the following term 

n random = $>i (Sf - 1/2) (0.3) 

i 

This term describes the interaction of the local spins with random magnetic fields. 

We now turn to the numerical calculation of the ground state energy as a function of 
the total magnetic flux E(&) which is the first step of our work. The ground state energy 
is obtained from a standard Lanczos algorithm ||. First of all, let us briefly describe some 
technical aspects of the method. The calculation is limited to relatively small system sizes 
N, since (i) the Hilbert space dimension grows exponentially fast with N, (ii) we have to 
average over many realizations of the disorder, (because of statistical fluctuations) and (iii) 
the disorder breaks the translation symmetry. The system sizes we were interested in, vary 

o 



from 6 to 20 sites and we have chosen to consider two different cases: half and quarter 
fillings. Let us remind that the Lanczos method consists in the construction of a tridiagonal 
matrix by applying iteratively the Hamiltonian on an initial random vector. By this way 
a basis of normalized vectors is defined as well as a set of values e n and b n given by 
the relationship, H^ n = 6 n _i^ n _i + e n \J/„ + b n+ i^ n+ i. Hence we construct by iteration a 
tridiagonal Hamiltonian matrix expressed in the \l/ n basis that we diagonalize to obtain the 
spectrum of the eigenvalues. This kind of process is rapidly converging. 

Let us, now, consider the calculation of the persistent current in such rings threaded by 
a total flux $. As usual the current is defined by 

where E(&) is calculated by exact diagonalization of the Hamiltonian. As well known, the 
flux can be gauged out from the Hamiltonian so that the presence of an Aharonov-Bohm 
flux through a ring is analogous to a twist in the boundary conditions ty(x + N) = ^(x)e 2m ^ . 
The spectrum and the persistent current have the flux periodicity of one. 

In figs. 1,2, /($) is plotted versus $ for a 16 site ring at electron density < n >= 0.5 
and < n >= 0.25 in (a) and (b) respectively, for different values of the disorder and the 
interaction. Fig.l corresponds to the 'ordered' interacting case W = and disorder W = 0.5 
is introduced in fig. 2. On fig.l, a discontinuity of /($) appears at $ = 0, for zero disorder 
(W = 0). Indeed, in the absence of disorder, translation symmetry is preserved and total 
momentum is a good quantum number. As a function of $, a crossing occurs between two 
lowest energy levels with different momenta. This crossing occurs at $ = for an even 
number of electrons or at $ = 0.5 for an odd number. This can be easily understood in 
the non interacting case where E(&) = — X^ ngf ^ cos(k n + 2irQ/N) and the subset of 
the electron momenta k n = 2irn/N is chosen in order to minimize the total energy. Since 
translation invariance is preserved in the presence of interactions, the discontinuity still exists 
for finite V. When disorder is introduced (W ^ 0) in fig. 2, the scattering potential lifts the 
degeneracy at the crossing point and hence leads to a continuous variation of the current. 
In fig. la (half-filling) we clearly observe the effect of the Mott transition on the currents: we 
notice that for V < 1, I is slowly varying with V, but when V > 1 a drop of the current 



appears. However, away from half- filling (see fig. lb), the current is not affected for moderate 
interactions. At W = and away from half-filling, the system is always metallic. 

In figs. 2a and b, the influence of the repulsive interaction is shown for a fixed impurity 
potential of magnitude W = 0.5 and for the same parameters as in figs.l (N=16, < n >= 0.5 
and < n >= 0.25). Clearly, in the half-filled case (fig. 2a), the repulsion tends to suppress 
the current even further. This is reminicent of the Mott localization which occurs in the 
pure system. With increasing W, /($) decreases as expected due to a stronger localization 
by the impurity potential. Such an effect has also been found in a ring of spinless fermions 



with long-range interactions 11 



More interesting is the effect of the interaction away from half-filling where no localization 
is expected in the absence of disorder. As seen previously in fig. lb the effect of the interaction 
for W = is extremely weak, because of the absence of Umklapp scattering. However, it is 
clear from fig. 2b that the interaction is much more effective in the presence of the disorder 
(W 0); while on fig. lb V had almost no effect, in fig.2b, for a relatively weak disorder, V 
leads to a significant decrease of the current. Such a striking influence of the interaction V 
is also seen in the metallic regime at half- filling (compare e.g. V = and 1 in figs, la and 
2a). 

At this point, we would like to describe more qualitatively the transition from the localized 
regime (insulator) to the ballistic one (perfect metal). As stressed by Scalapino and al. ||12|| , 
the Drude weight 7rD is a relevant parameter to characterize both of them. As originally 



noted by Kohn [13], the Drude weight nD can be calculated from the dependence of the 



ground state energy versus $ 



N ,d 2 E($), . . 

D = ^"fl^W (0 ' 5) 



where $ m = or 1/2 is the location of the minimum of E(<&). As mentionned earlier, $ m 
depends on the parity of the number of electrons. For an even electron number parity, we 
take D as the second derivative at $ m = 1/2. 

Note that for free electrons D =< n > /m where < n > is the density of the mobile 
charge carriers and m — l/2t is their mass. Generally D is given by |L4] : 



where m* is the effective mass of the carriers renormalized by the interaction. By this way, 
the D parameter can determine the different regimes. A perfect metal (ballistic regime) will 
be caracterized by a finite value of D. This corresponds to a persistent current I scaling as 
1/N. In the insulator, D vanishes exponentially as the size of the system goes to infinity, 
D oc e~ N ^, where £ is the localization length. As a check of the numerical calculations 
on the drude peak, we observed the correct finite size scaling for W = 0, at half-filling 
D ~ Dum + a/N 2 for V < 1 and D oc for V > 1 when N > 

Since we consider disordered systems, we have to average over many realizations of the 
disorder. The number of configurations we averaged over vary from 50 to 250, depending on 
the size of the Hilbert space and the filling. In fig. 3 (< n > = 0.5) and fig.4 (< n >= 0.25) D 
is plotted versus 1/N (N is the size of the ring). We first consider the half-filled case (fig.3). 
When V=0 and W=0, D goes to a finite value in the thermodynamic limit (Aim = I/ 71 ")- 
As long as W=0 (i.e. without disorder), the Drude weight is weakly affected by a small 
interaction V < 1, signature that the system remains metallic. But for V > 1 (here V=2) 
D decreases faster with N, the system becomes an insulator. We now turn to the effect of 
disorder. Once W ^ 0, we observe a tendancy towards localization for all V. With disorder 
in the system, the effect of the interaction is also to increase the degree of localization. In 
the light of the results given by fig.3 and 4, it appears that the role played by the interaction 
in presence of disorder (W ^ 0) is clearly different for < n >= 0.5 and for < n >= 0.25. 
On one hand, at < n >= 0.5 we see that the interaction V leads to a real decrease. For 
example, if we consider in fig.3 the case W = 2 and compare the data for V = and 2 we 
notice that the presence of the interaction reduces the localization length £ by more than a 
factor 3. On the other hand, at < n >= 0.25 (and fixed W) the value of D is less affected by 
the interaction V, as also observed on figs. 2a and 2b for the current. However, we note that 
the effect of V is significantly larger for W ^ than in the pure case W = 0. In summary, 
we do not observe any increase of D due to the competition between the interaction and the 
disorder. 

We finish this paper by a few remarks on the conductivity spectrum. The optical con- 



ductivity is given by , 

a(uj) = ttD5(uj) + a reg (oj) (0.7) 
where a reg (uj) given by the Kubo formula 

and j is the current density operator 

j = -it/2 J2i4 c i+ i exp{2i7r<S>/N) - he) (0.9) 

i 

All quantities in (9) are calculated at $ = $ m and E n are the excited manybody energies. 
The amplitude D of the Drude peak was calculated previously. We have explicitely 
checked the sum rule, 

/ a(u))duj = -(tv/2N) < 0\H kin \0 > (0.10) 
Jo 

where < 0\H kin \0 > is the groundstate expectation value of the kinetic energy. 

In fig. 5, (?{uj) is plotted versus uo for a 16 site ring at half-filling, with V = 0.5 and for 
W = or W = 1. In presence of the disorder we had averaged over 100 configurations. 
We clearly see, as expected, that in the metallic case W = 0, the contribution at non zero 
frequency is negligeable (ttD ~ 0.98 and only less than 1% of the weight is left at finite 
frequencies). However, when we introduce disorder, a strong absorption appears at non zero 
frequency with a peak around to = OAt, while weight is removed from uj = (ttD ~ 0.52). 

This broad absorption can be interpreted physically in the following way: if we assume 
that disorder localizes the wavefunctions in small finite size regions with a broad distribu- 
tion of volumes, this leads to a correspondingly broad distribution of finite characteristic 
frequencies. In other words, the localized electrons can oscillate in disconnected regions of 
different sizes. 

Let us summarize the main result of this paper. In our model for interacting electrons, we 
never observe an increase of the persistent current when interaction is switched on. At half- 
filling, the interaction induces a metal-insulator transition . The current strongly decreases, 
in qualitative agreement with the result of ref . |~l| • Away from half-filling, the effect of the 



interaction is much weaker in the absence of disorder. However when impurity scattering 
exists, the interaction plays again a crucial role and leads to an additional decrease of the 
current. This is because it is more difficult to move correlated electrons in a random potential 
than independent electrons. Our results are in discrepancy with those of J7|. It will be of 
interest to know if the qualitative results obtained in this paper still apply for a multichannel 
ring. 
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FIGURE CAPTIONS 

Figure 1 

Current /($) versus $ for a 16 site ring at < n > = 0.5 (fig.la) or at < n >= 0.25 (fig.lb), 
fixed V = and W = 0, 0.5 and 1. 
Figure 2 

Current /($) versus $ for a 16 site ring at < n >= 0.5 (fig.2a) or < n >= 0.25 (fig. 2b) ,at 
fixed W = and V = 0, 1 and 2. 
Figure 3 

Scaling of D at half-filling (< n >= 0.5). D vs 1/N for different V and W. 
Figure 4 

Scaling of D at quarter filling (< n >= 0.25). 
Figure 5 

Total conductivity a(u) for a 16 site ring at < n >= 0.5 and (a) with or (b) without disorder 
(W = and W = 1) for a fixed V = 0.5. For W = tcD ~ 0.98. For W = 1, 7rD raer . ~ .52 
and we averaged over 100 realisations of the disorder. 
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